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RESEARCH MEMORANDUM 

CALCULATIVE METHOD FOR ESTIMATING THE INTERFERENCE PRESSURE 
FIELD AT ZERO LIFT ON A SYMMETRICAL SWEPT-BACK WING 
MOUNTED ON A CIRCULAR CYLINDRICAL BODY 
By Jack N. Nielsen and Frederick H. Matteson 


SUMMARY 

An approximate method is presented for calculating the inter- 
ference pressure distribution at zero lift of a symmetrical swept— 
back wing mounted on a circular cylindrical body with the chord plane 
of the wing passing through the axis of the body. The method 
employed is based on linear theory and is an approximate first-order 
method. An example illustrating the method is presented for a wing- 
body combination consisting of a circular cylindrical body and an 
untapered wing swept back 60°. The wing considered in the example 
has a double-wedge section with the maximum thickness at the 
midchord, a chord equal to twice the body diameter, and a spanwise 
distance from the wing-body juncture to the wing tip of two body 
diameters. The effects of a body nose on the wing pressure distribu- 
tion are not considered in determining the pressure drag of the wing 
in combination. It is found that the pressure drag of the wing in 
combination is less than that of the two exposed half— wings joined 
together . 


INTRODUCTION 

At subsonic speeds, wing— body interference has important effects 
on the aerodynamic characteristics of an airplane. At supersonic 
speeds, wing-body interference assumes even greater importance because 
of the relatively large bodies and small wings usually specified for 
aircraft designed for flight at such speeds. Swept>-back wings 
probably have greater interference effects at supersonic speeds than 
unswept or swept— forward wings because greater areas of swept— back 
wings lie within the region of influence of the wing— body juncture. 

Past work on wing-body interference at supersonic speeds has been 
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largely confined to lift and moment interference. Using the concepts 
developed by R, T. Jones in his low-aspect— ratio, triangular— wing y 

theory (reference 1), Spreiter (reference 2) has determined the 
effect of a pointed, slender "body of revolution on the lift— curve 
slope and center-of— pressure location of a wing-body configuration. 

Browne, Friedman, and Hodes (reference 3) , using the methods of 
conical supersonic flow, have determined the loading on a triangular- 
wing conical— body combination with a common apex. Ferrari (refer- 
ence 4) has analyzed the lift and moment interference between a flat 
rectangular lifting surface and a pointed body of revolution. 

There has been little work on the effects of wing thickness on 
the interference pressure distributions of wing-^body combinations at 
supersonic speeds. It is the purpose of this report to present an 
approximate calculative method for determining the pressure field due 
to interaction between a wing with a double-wedge section that is 
swept behind the Mach cone and the circular cylindrical portion of 
the body of revolution in which it is mounted. The combination is 
analyzed at zero angle of attack with the chord plane of the wing 
passing through the body axis. The interferehce pressure distribu- 
tions, acting at several spanwise positions of the wing in cambina- i 

tion, are given in the report together with the wing interference 
pressure drag. 

y 

The emphasis in the analysis is placed on swept— back wings, the 
surfaces of which are composed of planes. The entire analysis, based 
on the linearized theory of supersonic flow, is for a free-stream 
Mach number of */2. The method for applying the results of the 
analysis to other Mach numbers is given at the end of the text. The 
approximate nature of the solution arises from the mathematical com- 
plexity of the complete interference problem. One of the components i 
of the interference pressure field requires the solution of an 
integral equation. Since the solution was not found, it was neces— , 

sary to resort to a tedious numerical summation process. 


SYMBOLS 

c wing chord 

Ci, C 2 , C 3 factors for determining wing pressure coefficients 

ca section drag coefficient 


d 


body diameter 
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k 1# kg, k 3 

m 

Mo 

n i> n 2j b- 3 , n* 
P 

5o 

P 


factors for determining average normal velocity 
at control area due to set of rectangular wedges 

cotangent of sweep angle of line pressure source 

free— stream Mach, number 

averages used in determining section drag coeffi- 
cient 

local static pressure 
free— stream static pressure 

pressure coefficient ft?) 


So 


r.p. 


s 


so 

t 

u,v,w 

Uw* Vw, Wv 




free— stream dynamic pressure 
real part of a complex function 

semispan of wing measured from wing— body juncture 
to wing tip 

semispan of rectangular wedge 

maximum thickness of double-wedge section 

perturbation velocities of line pressure source 

perturbation velocities for rectangular wedge in 
x w> Yv) z w directions, respectively 

finite velocity normal to body (positive outward) 
due to leading— edge sources after reduction of 
infinite sidewash velocities to finite velocities 


vn average finite normal velocity for control area 

v n ’ average velocity normal to control area due to 

rectangular wedge 


V 0 


free— stream velocity 
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x»,y',z« 
x w> yv, z w 


(S) s 

'dz\ / dz\ /dz\ 

^VaWbWc 


fcX-I 

(s> 



longitudinal, lateral, and normal coordinates for 
line pressure source or sink with source or sink 
in z=0 plane and the origin of the source or 
sink at x=0, y=0 

(The positive x direction is downstream, the 
positive y direction is from the x axis toward 
the source or sink, and the positive z direction 
is upward.) 

oblique coordinates (x’=x-my, y’ =y— mx, z’ =z //i-m 2 ) 

longitudinal, lateral, and normal coordinates for 
rectangular wedge with chord plane of wedge 
coinciding with Zv=0 plane 
(The positive x v direction is downstream, the 
the positive y v direction to the right for an 
observer looking upstream, and the positive 
z w direction away from the body.) 


slope of streamwise section due to line pressure 
source or sink 

(Positive for a source and negative for a sink. ) 

V 


streamwise slope of rectangular wedges correspond- 
ing to strips A, B, and C, respectively 


streamwise slope of rectangular wedge correspond- 
ing to strip A and the region of influence which 
starts at the forward edge of control area A— I 


slope of rectangular wedge in streamwise direction 


ratio of average finite normal velocity to free— 
stream velocity for control areas of strips A, 
B, and C, respectively 


change in ratio of the average finite normal 
velocity to the free— stream velocity due to a 
given ring of rectangular wedges 
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ratio of average finite normal velocity to free— 
stream, velocity for control area A— I 


cp potential of the perturbation velocities of a line 

pressure source or sink 

9 polar angle in planes perpendicular to "body with 

origin at intersection of plane with body axis 
(The 0=0 plane corresponds to the right half- 
wing for an observer facing the oncoming stream 
and 0 is positive counterclockwise.) 

it r\, t, longitudinal, lateral, and normal coord ina tes with 

the origin at the leading edge of the wing— body 
juncture and the £=0 plane coincident with the 
chord plane of the wing 

(The positive £ direction is downstream, and the 
positive £ direction is up. The positive Tj direc- 
tion is taken to the right of an observer facing 
the oncoming stream for the right half— wing and is 
taken to the left for the left half— wing. ) 


ANALYSIS 

Components of the Interference Pressure Drag 

The interference pressure drag for a wing-body combination is 
the difference between the drag of the wing-body combination and the 
sum of the drags of the body alone and the wing alone. The wing alone 
is used in this report to mean the two exposed half— wings joined 
together. The interference pressure drag arises from the following 
pressure distributions: (1) The pressure distribution an the body in 
combination minus that on the body alone, and (2) the pressure distri- 
bution on the wing in combination minus that on the wing alone. This 
report considers a symmetrical wing mounted on the cylindrical portion 
of a pointed body of revolution. The longitudinal axes of the body 
and wing ard in the stream direction. Since the sides of the body in 
the region influenced by the wing are then parallel to the flow, the 
body has no interference pressure drag. Thus, the interference 
pressure drag is the difference in the drags of the wing in comb ina — 
tion and of the wing alone. The methods of references 5 or 6 can be 
used to determine the drag of the wing alone, but no method has 
hitherto existed for determining the pressure distribution and drag 
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of the wing in combination. It is the purpose of this report to 
present such a method. 


Components of Pressure field on Half— Wing in Combination 

In the following a n alysis of the pressure distribution and drag 
of either half— wing in combination, the pressure field acting on the 
half— wing is subdivided into four components: (1) That due to the 

part of the body forward of the region of influence of the wing, 

(2) that due to the given half— wing, (3) that due to the opposite 
half— wing, and (4). that due to interaction between the half— wings 
and the part of the body within their region of influence. This 
subdivision of the pressure field acting on the half— wing in combina- 
tion is convenient for the purposes of the ana lysis. However, once 
the method of the analysis has been established, a slightly different 
point of view will be adopted. 

Consider now methods for determining the various components. 

The first component can be determined by any of the methods available 
for determining the pressure field of bodies of revolution, Buch as 
the method of characteristics or the method of reference 7. In 
general, this component is different for each body shape and no 
further account will be taken of it in the analysis. The second and 
third components can be determined by the method of reference 5. 

This method is considered briefly in this report since It is the 
starting point of the analysis to obtain the fourth component, which 
is the one to which the analysis is principally devoted. It Bhould 
be noted that the part of the body forward of the region of Influence 
of the wings will cause only small perturbation velocities at the wing 
and will have only a second-order effect on the fourth component of 
the pressure field. 


Pressure Field on Half— Wing Due to Itself 
and Opposite Half— Wing 

The line pressure source solution of the linearized equation of 
supersonic flow is fundamental in determining the pressure field on 
either half— wing due to itself or to the opposite half— wing. Jones 
(reference 5) giveB for the streamwise perturbation velocity u of 
a line pressure source swept behind the Mach cone 
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The. velocity potential for the lin e pressure source is obtained by 
integration of the foregoing equation along a path parallel to the 
x axis. This gives 


<p = f udx = — f—') -fy* cosh -1 f —=£== - 

J_co ic \dx/ s l \ ,/y'2+z'2/ 


y a/ 1 -^d 2 oosh~i 


fe) 


+ Z* COS -1 


yy 1 +z g ' 1 \ 

,/(y2+z2)(y>2 + z«2) J J 

( 2 ) 


The lateral and vertical perturbation velocity components are 
obtained by differentiating equation (2) as follows: 

V = |2 = Xfl ($£) [— cosh- 1 ( 7 -fl - cosh- 1 ( -r*—') 

by it \dx/ B L *A-^ 2 V Vy ,2 +z ,2 y V-Zp+z 2 / 

(3) 


w = |2 « 2a /jto'N cos- 
* ^dxy/s 


dz 


yy 1 +z 2 


L V ( y^z 2 ) ( y 1 s +z * 2 ) J 




To illustrate the nature of the downwash and sidewash, these 
fields have been calculated from equations (3) and (4) for a sweep- 
back angle of 60° and a free— stream Mach number of ,/5. The res ul ts 
of the calculations are presented in figures 1 and 2. These figures 
show the downwash and sidewash patterns in any plane perp endi c ular 
to the axis of the Mach cone of the source. Since the flow is 
conical, the downwash and sidewash velocities depend o nly on y/x 
and z/x, the ratios which determine the direction of any ray from 
the apex of the Mach cone. The line source starts at the apex of 
the Mach cone and pierces any plane perpendicular to the axis of the 
Mach cone at a point y/x=0.577 and z/x=0. Figure 1 reveals that 
upwash exists above the plane of the source, and downwash exists 
beneath it. Between the line source and the axis of the Mach cone, 
the upwash and downwash velocities are constant just above and below 
the zero pl an e and their magnitudes are equal. Thus, behind the 
line source there exists, in effect, a symmetrical wing with a wedge- 
shaped streamwise section. 
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(In the linearized theory of supersonic flow, the surfaces of the 
wing are displaced only infinitesimally from, the z=0 plane.) In 
the rest of the z=0 plane, the downvash velocity is zero. The 
sidevash pattern of figure 2 is symmetrical above and below the z=0 
plane. There is positive infinite sidevash velocity along the line 
pressure source and negative ■ infinite sidevash velocity along the 
axis of the Mach cone. 

The fact that a line pressure source or sink produces zero 
dovnwash velocity everywhere in the zero plane except .directly behind 
it simplifies the simulation of complicated wings by the superposition 
of line sources and sinks. In fact, the pressure field of any wing 
which is symmetrical above and below the z=0 plane and the surfaces 
of which are composed of planes can easily be simulated by a finite 
number of line sources and sinks. For instance, the pressure field 
of an untapered, swept— back half— wing of symmetrical double— wedge 
section with the maximum thickness at the midchord would be formed 
by the source— sink system shown in figure 3* The leading-edge and 
trailing— edge sources are of equal strength, but the midchord— line 
sink is twice as strong as each of the sources. A source and two 
sinks are introduced to form a tip. In determining the pressure 
field of the half— wing from its source— sink system, the u perturba- 
tion velocity for each source or sink is determined from equation (l) 
for a given point on the wing, taking into account the fact that the 
starting point of each source or sink is taken as the origin of the 
x,y,z coordinate system in determining its particular value of u. 

The pressure coefficients of the line sources and sinks are then 
determined from the following equation which is based on linear theory 


p -- 2 0t) (5) 

Finally, the pressure coefficients are summed to give the pressure 
coefficient of the half— wing at the particular point. 

In deter mining the drag of the half— wing from the pressure field, 
it is not necessary to consider the contributions of the tip source- 
sink system to the wing pressures when the Mach lines of the tip 
source and sinks on a given half-wing intersect the trailing edge of 
this half— wing. In this instance, which is frequently the case, the 
tip source— sink system changes only the distribution of drag within 
its region of influence but does not change the over-all drag. (See 
references 8 and 9») 

The deter mina tion of the pressure field on one half— wing due to 
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the other half— wing is simplified "by the fact that one half— wing does 
not produce any downwash in the plane of the other half— wing (fig. 1). 
Hence, the pressure field due to the one can "be considered to act 
directly on the other. Figure 4 presents an example of the inter- 
action of the pressure fields of two half— wings separated hy half a 
chord length. For the configuration shown, the leading— and trailing- 
edge sources, as well as the midchord sinks, all contribute to the 
pressures on the opposite half— wing. In computing the contributions, 
it should be noted that the value of y to be used in determining y* 
in equation (1) will be negative. In this particular case, the tip 
system of one half— wing does not contribute to the pressures on the 
opposite half— wing. 


Pressure Field on Half-Wing Due to Interaction 
Between Half— Wings and Body 

The source— s ink system which represents the half— wings will 
induce both downwash and sidewash against that portion of the body 
within their region of influence. However, the velocity normal to 
the body must be zero everywhere. This constraining action of the 
body affects the wing pressures and gives rise to the fourth component 
of the pressure field acting on the half— wings in combination. To 
determine this component of the pressure field, it is necessary to find 
solutions of the linearized equation for supersonic flow which will 
cancel the normal velocity induced at the body surface by the half- 
wings and which will have no contribution to the downwash in the plane 
of the wing. It is apparent that a method for determining that part 
of the fourth component associated with the leading— edge sources will 
be equally applicable to determining those parts due to the midchord— 
line sinks and the trailing— edge sources. For certain restricted con- 
ditions, such as for wings of very low aspect ratio or for very low 
supersonic Mach numbers, the tip source-sink systems may have an 
effect on the fourth component of the pressure distribution on the 
half— wing. However, the possible effect of the tip systems is ignored 
in the rest of the analysis, and only the leading-edge sources are 
considered. 

Reference to figures 1 and 2 shows that infinite normal velocity 
can arise from the sidewash only. Infinite sidewash velocity is 
induced against the body along the intersections of the body with the 
wing-chord plane. Finite downwash and sidewash velocities are 
induced everywhere else on the body. It Is convenient to reduce the 
infinities in the normal velocity field to finite values and then to 
cancel the finite residue. As will be shown next, this infinite 
sidewash can be reduced by a source— sink system. 
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Reduction of Infinite values of sidewash velocity to finite 
values by source-sink syst em.- The infinite sidewash velocities 
normal to the body in the £ =0 plane due to the leading— edge sources 
can be reduced to finite velocities by means of the source— sink 
system shown in figure 5- Image source II cancels identically the 
infinite sidewash velocities due to source I, and image source V 
cancels identically those due to source IV. The image sources, how- 
ever, produce infinite sidewash velocities along themselves (fig. 2) 
and produce downwash velocities behind themselves which in terfere 
with the wing boundary conditions. These two conditions require the 
addition of image sinks III and VI, the origins of which lie on the 
image sources inside the body. As shown in figure 5, the origins of 
the image sinks are displaced from the origins of the corresponding 
sources. (Otherwise the sinks would identically cancel the sources.) 
It is noteworthy that the infinite sidewash velocities along the 
image sources are not identically canceled by the infinite sidewash 
velocities of the corresponding image sinks because of the displace- 
ment of their origins. The resultant sidewash velocity at the body 
is, however, finite. To accomplish the reduction of the infinities 
in a satisfactory manner, the origins of the image sinks must be 
close to the origins of the .corresponding image sources. This point 
will be considered in detail in the subsequent illustrative example. 

The contribution to the pressure distribution of the wing due 
to reduction of the infinite sidewash is readily obtained. It is the 
sum of the pressure coefficients due to the sources II and V, and 
sinks III and VI as determined from equations (1) and (5). Again, 
the origin of the x,y,z coordinate system is taken at the starting 
point of the particular source or sink for which the calculation of 
the pressure coefficient is being made. 

After the infinite sidewash velocities have all been reduced to 
finite values, there remains a distribution of finite velocity normal 
to the surface of the body. The part ■ of the fourth component of the 
wing pressure field arising from canceling these finite normal veloc- 
ities will be only approx ima tely determined. 

Approximate cancellation of the finite normal velocities by 
rectangular wedges .— The finite velocity normal to the "body at any 
point after the cancellation of the infinite sidewash may be deter- 
mined with the aid of equations (3) and (4). From these equations 
the sidewash and downwash due to sources I, II, IV, and V and sinks 
III and VI are determined, and are then summed by the following equa- 
tion to obtain the finite velocity normal to the body: 
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i=VI i=VI 

?5- I (t)i COB9+ I (jk) ± Bl * e (6) 
1=1 1=1 


The right-hand sideline of the body for an observer facing up- 
stream. corresponds to 0=0°, and the left-hand sideline to 0=l8O°. 

A positive v n /V o denotes flow out of the body with v positive to 
the right (looking forward) and w positive upward. Only the normal 
velocity distribution for values of 0 from 0° to 90° need be con- 
sidered since the flow is similar in every quadrant of the body. 

The surface of the circular cylindrical body is now divided into 
control areas , the number of ■which determines the accuracy of the 
solution. These control areas are formed by the intersections of 
longitudinal strips with a number of rings as shown in figure 6. 

Eight longitudinal strips and eighteen rings were selected in the 
present analysis to provide a reasonable compromise between accuracy 
and work required. Strip A (fig. 6) has a perpendicular tangent to 
the plane of the wing, and the value of 0 for this strip varies 
from -22.5° to +22.5°. The rings are 0.2d wide and start at g=0, 
where the leading-edge sources intersect the body. 

The next step in the analysis is to present a method for cancel- 
ing the average flow induced through each control area by leading-edge 
sources. The finite normal velocities at every point on the body 
surface. can, in principle, be identically canceled by a distribution 
on the surface of point sources of the type used by Puckett in refer- 
ence 6. For a plane distribution of these sources, the local source 
strength is proportional to the local induced vertical velocity. 
However, when the sources are placed on a curved surface, this simple 
relationship no longer holds, and the determination of the surface 
distribution of the source strength to cancel an arbitrary surface 
distribution of normal velocity involves the solution of an integral 
equation. Since the solution of this equation was not obtained, it 
was necessary to use an approx ima te method. 

The approximate method used to cancel the average finite normal 
velocity at each control area employs the solution for a rectangular 
v'edge of infinite chord. The rectangular wedge is formed by the 
source-sink system shown in figure 7* The sources and sinks all have 
no sweep and correspond to m=oo. For m=eo, the u, v, and w velocity 
components of a source given by equations (l), ( 3 ), and (^) reduce to 
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The u, y, and w components for the rectangular wedge with the 
origin as shown in figure 7 are 



The rectangular wedges are placed on the sides of a regular 
octagonal prism inside the body as shown by figure 8. From symmetry 
considerations, it is apparent that the wedge angles, or slopes, of 
the rectangular wedges to be used are identical for the two vertical 
sides of the prism, the two horizontal sides, and the four inclined 
sides. It is also apparent that no flow will be induced across the 
plane of the wing, and therefore the wing boundary conditions are 
satisfied. The rectangular wedges would induce infinite normal 
velocities at the cylindrical surface if the octagonal prism were 
inscribed in the circular body. To avoid this difficulty, the 
octagonal prism has arbitrarily been inscribed inside a cylinder 


v 
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concentric with, the body and having a diameter 0 . 966 .. Rectangular 
wedges are quite suitable for the present purpose since it is neces- 
sary only to satisfy the boundary conditions at the body surface and 
not at positions inside the body. 

A step-by— step process is used in determining the slopes of the 
rectangular wedges. First, rectangular wedges are placed on the 
octagonal prism so that their region of influence on the body begins 
at the 5=0 plane (fig. 6). From equations (ll) and (12) for the 
dowhvash and sidewash, the ratio of average normal velocity to free— 
stream velocity vn’ /Vo through any control area due to the rectan- 
gular wedge on side A of the octagonal prism is then determined. This 
average normal velocity ratio is dependent on a factor k which is 
different for each control area as follows: 



where (dz/cbc)A refers to the slope of the rectangular wedge corre- 
sponding to strip A. The values of k for all control areas above 
the horizontal plane of symmetry have been determined. The value of 
the normal velocity was calculated with the help of equations (ll) 
and (12) for a number of points on the body and yn.’/v'o for each 
control area was determined by a process of cross— plotting and 
integrating with a planimeter. Finally, the values of k were 
determined so that equation (13) 'was fulfilled. Table I gives the 
factor k for the various control areas due to the rectangular 
wedge on side A of the • octagonal prism. The table presents results 
for only half the body since the results are similar above and below 
the plane of the wing. It is apparent that the factors in table I 
may be used to deter min e the average normal velocity through any 
control area due to any rectangular wedge. 

With the aid of the factors in table I, the slopes of the 
rectangular wedges required to cancel the average finite normal 
velocity ratio v n /Vo in every control area of ring I can be deter- 
mined. Let (dz/dx)A_x be the slope of the wedges which lie on sides 
A and E of the prism (fig. 8) and the regions of influence of which 
start on the body at 5=0 (fig. 6). Let (dz/dx)B_x ^e corre- 
sponding slope of the wedges for sides B, D, F, and H and (dz/dx)o_i 
the slope of the wedges for sides C and G. Now let ( v n /Vo)^_j be 
the ratio of the average finite normal velocity to the free— stream 
velocity for control areas A— I and E— I , (vn/v 0 )p_x that for control 
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areas B— I, D-I, E-I, and H-I, and (v n /V 0 )Q_j that for control areas 
C— I and G— I. Consider now control area A— I. On the "basis of the 
factors in table I, the ratio of average normal velocity to free— 
stream velocity induced through this control area by the wedge of 
side A corresponds to 0.8l3(<lz/dx)A_. Since a wedge of given slope on 
side A will induce the same normal velocities through control area 
B— I as a wedge of identical slope on side B will induce through 
control area A— I, it follows from table I that the ratio of average 
normal velocity to free— stream velocity induced through control area 
A— I by the wedge on side B is 0.069(dz/dx)B— I. By symmetry, the wedge 
on side H will also have t'he contribution 0.069(dz/dx)B_j* The wedges 
on sides C, D, E, E, and G have no contribution because the control 
area A— I is outside their regions of influence for the Mach number 
and size of control area selected in the analysis. Since the average 
velocity normal to control area A— I is to be zero 

+ 0.813 (—^ +2(0.069) =0 ( 14 ) 

A—I . 


Similarly, for control areas B— I and C— I, there is obtained 


('"M + 0.813 ffe') + 0.069 + (#*) 1-0 

\Vo/b-i V^x/e-i L a _i c-i -1 (15) 

and 

( — S ) + 0.813 + 2(0.069) ( 4^ ) » o (16) 

\Vowc— I \dx/Q_ 1 \dxs 3—1 


Ehe simultaneous solution of equations (14) to (16) gives the slopes 
of the rectangular wedges, 


\ Y ° /B-I 


(£;« - - + °- 2i K - °- oi8s5 ( 

(s)« ■ °- 1075 [(*)« + (^L ]- i - ser (^) I 


(18) 
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= - 0.01825 


In 

Vo 



+ 


0.2151 




i. 2W (a) 


C-I 


(19) 


These values of (dz/dx)^_i, (dz/dx)s_i, and (dz/dx)c— I are such that 
the average normal velocity for each control area of ring I is zero. 

I 

Once the slope of the rectangular -wedges required to cancel the 
average finite normal velocities in the control areas of ring I are 
determined from equations ( 17 ) to ( 19 ) , the values of v^/Vq for all 
the control areas "behind ring I must he adjusted for the contributions 
due to these rectangular wedges. The adjustment need be made only 
for' strips A, B, and C because of symmetry. Consider, for example, 
control area B— VII. Due to the wedge on side A of the octagonal 
prism, the ratio of average normal velocity to free— stream velocity 
for this control area is 0.2l8(dz/dx)A— l> using the factors of 
table I. Consider now the contribution of the wedge on side E. From 
symmetry, the effect on control area B-VII of the wedge on side E is 
the same as the effect on control area D— VII of the wedge on side A. 
Thus the contribution of the wedge on side E to the ratio of average 
normal velocity to free-stream velocity for control area B— VII iB 
0 . 165 ( dz /dx)^—! . The total contribution of sides A and E is thuB 
(0.218 + O.I 65 ) (dz/dx)A— I- In a similar manner, using the factors 
of table I, it can be shown that the total contribution for the wedges 
of sides B, D, E, and H is (O .835 + 0.153 + 0.180 + 0.153) (dz/dx);B_i 
and that the sum contribution of sides C and G is (0.218 + O.I 65 ) 
(dz/dx)c— 1« From this example, it can be seen that the contribution 
to the ratio of average normal velocity to free-stream velocity due 
to the rectangular wedges of a given ring can be represented as 


A 



) 



( 20 ) 


The values of k x , k 2 , and k^ have been determined in the foregoing 
manner for all the control areas of strips A, B, and C. The factors 
k x , k 2 , and h 3 , due to the rectangular wedges of ring I, are 
presented in table II . When these factors are used, the values of 
Vn/Vo for the control areas behind ring I are readily adjusted for 
the normal velocity induced by the rectangular wedges of ring I. 

The foregoing process Is repeated in the stejn-by— step solution. 
A second set of rectangular wedges is placed on the octagonal prism 
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so that the region of influence on the body begins at the bo undar y 
between rings I and II at £=0.2d. The average finite normal veloc- 
ities for the control areas of ring II are then canceled in the same 
manner as for ring I, using equations (IT), (18), and (19), but sub- 
stituting the quantities (v n /v 0 ) A _j J} (v n /V'o) B _ II , and (v n /V 0 ) c _ II 
for (yq/Vo)a-I 0WVo);b-I> and (vn/Vo)c_i* respectively, and solving 
for the quantities (dz/dx)^— II, (dz/dx)s_jx, and (dz/dx)c_n. The 
values of v n /Vo for all control areas behind ring II We adjusted 
as previously, using equation (20) and the factors of table II. In 
using the factors of table II, however, it must be borne in mind that 
ring II is farther downstream than ring I. Thus, for instance, the 
average normal velocity induced at control B-VII by the wedges of 
ring II would be the same as that induced at control area B-VI by the 
wedges of ring I. Thus, for the control area B-VII, using the factors 
in table II, 


A 


/vin\ 

vW B-VII 


0.433 


A— II 


1.4-12 


(£)*.* 


+ 



In the foregoing manner, the slopes of the wedges for ring II are 
determined, and the average normal velocities for all control areas 
behind ring II are adjusted for the normal velocities induced by the 
wedges of ring II. 


The procedure used in determining the slopes of the wedges for 
ring I and adjusting the control areas behind ring I for the induced 
normal velocities of the wedges of ring I has been applied to ring II. 
It Is clear that this procedure must also be applied to each additional 
ring in the order III, IV, V, etc., through the last ring, that Is, 
through the last ring the region of influence of which intersects the 
wing. With regard to the procedure, it is desirable to emphasize the 
use of table II. Consider the contribution to the ratio (vn/Vo) at 
control area B— VII due to the wedges of ring VI. This contribution 
will be the same as that at control area B-II due to the rectangular 
wedges of ring I. Thus, for the effect of ring VI on control area 
B-VII 


A 


= 0.189 TfO + 0 . 88 j(- 

\V 0 / B— VII \&x/a-VI v 


A— VI 


&0 + 0.189 (V) 

^'B-VI \dx/ q_yj 

( 22 ) 
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The procedure to determine the slopes of the rectangular wedges 
is now summarized. Equations (17) , (18) , and (19) are solved to 
obtain the slopes of the wedges for ring I. The values of (vn/v o) 
for all control areas behind ring I are adjusted for the normal 
velocity induced by the wedges of ring I using equation ( 20 ) and the 
factors of table H. The slopes of the rectangles of ring II are 
then deter min ed using equations similar to equations ( 17 ), ( 18 ), anA 
( 19 ) with identical coefficients. The values of v n /Vo f° r a li 
control areas behind ring II are adjusted 'using equation (20) and the 
factors of table II, and the process is repeated through the last 
ring tha t can influence the wing pressures. 

From the slopes of the wedges, that part of the fourth component 
of the pressure field acting on the half-wing in combination which is 
due to canceling the average finite normal velocity at each control 
area can be calculated. This calculation completes the determi n ation 
of the pressure field acting on the half— wing in combination. The 
pressure coefficient at any point due to any one of the rectangular 
wedges can be calculated with the help of equations (5) and (10). It 
can be easily shown that the pressure coefficient at any point on the 
wing due to the first set of rectangular wedges can be expressed in 
the form 



The contributions of the other sets of rectangular wedges to the 
pressure at a given point can be similarly expressed and calculated. 

The contributions of all of the sets are summed at a number of points 
at a given spanwise station to yield the Bection pressure distribution. 
The details of calculating the pressure coefficients will be given in 
the example. 


ILLUSTRATIVE EXAMPLE 

’Tbp pressure field act ing on a half-wing in combination has been 
subdivided into four components in a particular ma nn er to facilitate 
the analysis. The fourth component, that due to interaction between 
ha lf— w ing s and body, was further subdivided into two parts. With the 
deter mina tion of its various components, the pressure distribution on 
a half-wing is fully determined. The part of this pressure distribur- 
tion due to interference will be considered in this illustrative 
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example. The over-all interference pressure distribution of the 
half— wing in combination is defined for the purpose of this report to 
he the sum of that component due to the opposite h a lf-w ing and that 
component due to interaction between the half— wings and body. The 
component due to the body nose is not considered for reasons already 
given, and the component due to the effect of the half-wing on itself 
is not an interference pressure distribution. The over— all inter- 
ference pressure distribution is considered in two parts to facili- 
tate the fairing of curves required in this illustrative example. 

The first part of the overall interference pressure distribution iB 
the sum of the pressure distribution due to the opposite half— wing 
and the pressure distribution due to reducing the infinite sidewash 
velocities to finite velocities. The second part is the pressure 
distribution due to canceling the average normal velocities at each 
control area by rectangular wedges. 

In this illustrative example the interference pressure distri- 
bution at zero lift and the interference pressure drag are calculated 
at M 0 = \/2. for a 60° swept— back wing mounted on a circular cylindri- 
cal body . The dimensions of the wing— body configuration considered 
are given in figure 9 . An untapered svept-back wing has been chosen 
which has a symmetrical doub le— wedge section with the maximum thick- 
ness at the midchord, and the body diameter has been taken equal to 
one— half the wing chord. The over-all interference pressure distri- 
butions due to the lead in g— edge sources have been determined for the 
stations shown in figure 9- The spanwise distribution of pressure 
drag corresponding to the over-all interference pressure distributions 
has been calculated. Finally, the interference pressure drag, as 
defined in the analysis, has been determined. 


Wing Pressure Field Due to Leading— Edge Source of 
Opposite Half— Wing and to Canceling Infinite 
Sidewash of Leading— Edge Sources 

The pressure distribution at the wing-body juncture of the 
right half— wing (for an observer looking forward) due to the leading- 
edge source of the left half— wing and due to reducing the infinite 
sidewash velocities of the leading-edge sources to finite velocities 
is shown in figure 10. Only the interference pressure distribution 
due to the leading-edge sources is shown for purposes of clarity. 

The corresponding pressure distribution for the midchord— line sinks 
is double that for the leading— edge sources, is negative, and is 
shifted a half of a chord rearward. The corresponding pressure dis- 
tribution for the trailing-edge sources is identical to that for the 
leading— edge sources but is shifted a chord length rearward. 
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In figure 10 the contribution of the lea ding -edge source of the 
left half—' wing is that due to source IV and the contribution due to 
reducing the infinite sidewash velocities of the leading— edge sources 
is represented by the sum of sources II and V plus sinks III and VI. 

The total effect is II + III + IV + V + VI as shown. 

These pressure distributions present several noteworthy features: 
First, the wing pressures due to reducing the infinite sidewash veloci- 
ties are all finite since the infinite pressure peak of source V is 
canceled, although not identically, by the infinite pressure peak of 
sink VI without forming a cusp. A cusp near the leading edge arises 
in the sum pressure distribution (II + HI + IV + V + VI) due to 
sink III. The position of this cusp depends only on the position of 
the origin of sink III, which is arbitrary. As the origin of sink III 
approaches that of source II, it is apparent that the effect of the 
.cusp is made smaller. Two cusps are introduced into the sum pressure 
distribution at the intersections of the Mach . l in es of source V and 
sink VI with the wing-body juncture. As the origin of sink VI 
approaches the origin of source V, the difference in the heights of 
the cusps decreases as well as the chordwise distance between them. 

For reasons that will subsequently be discussed, it is desirable to 
minimize the effect of the cusps in the final pressure distribution 
by locating the origins of sinks III and VI close to the origins of 
sources II and V, respectively. It was found that positions of the 
origins of sinks III and VI of either 0.05<1 or O.lOd from the wing— 
body juncture gave good results in the present example, as will be 
pointed out. If the origins of the sinks are located too close to the 
origins of their respective sources, however, the finite velocities 
normal to the body (v n ) would be too large along the wing-body junctures 
and would impair the accuracy of the method. 

The pressure distribution In figure 10 represents only part of the 
interference pressure distribution at the wing-body juncture due to the 
leading— edge sources . The second part arises from canceling the finite 
velocities normal to the body. As previously mentioned, a distribution 
of point sources of the type used by Puckett scattered over the body 
could, in principle, cancel the values of VnAo at every point on the 
body. The pressure distribution at the wing-body juncture resulting 
from this source distribution would have cusps that would directly 
offset those of figure 10. However, in the present approximate method 
of obtaining the second part of the interference pressure field, the 
offsetting cusps will not be obtained. 

The manner in which the pressure distribution, due to the leading- 
edge source of the opposite half— wing and to reduction of the infini te 
sidewash velocities of the leading— edge sources, varies in the spanwise 
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direction is shown in figure 11. Near the leading edge of the wing- 
body juncture (t)/s=0} the pressure coefficient is large. This is to 
be expected, as will be pointed out. The pressure coefficients are' 
all zero at the forward Mach line for all stations except the wing- 
body juncture, " and the first cusp reduces in severity as the spanwise 
distance increases. Although the double cusp near the midchord of 
the inboard station reduces in severity as the spanwise distance 
increases, it does not do so as rapidly as the cusp behind the 
leading— edge Mach line. The parts of the pressure distributions 
between the leading edge and the leading— edge Mach line have been 
included in the figure because these distributions are moved downstream 
in determining the corresponding contributions of the midchord— line 
sinks and the trailing-edge sources on the wing. 


Wing. Pressure Field Due to Approximate Cancellation 
of the Finite Velocities Normal to the Body 
Due to Leading-Edge Sources 

The part of the interference pressure field due to approximate 
cancellation of the finite velocities normal to the body has been 
determined only for the leading— edge sources as before. First, the 
ratio of the average finite normal velocity to free— stream velocity 
v n /Vo for each control area after cancellation of the infinite Bide- 
wash was determined. The values of v/Vo ‘and w/v 0 for sources I, 
II, IV, and V and sinks III and VI (fig. 5) were calculated for 
values of 9 0°, 10 °, 30°, 60°, and 90°, and for distances |/d 

of 0, 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0, using equations (3) and ( 4 ). 
The values of v n /v 0 were then calculated from equation (6) . The 
values of yn/Vo were plotted against |/d for each value of 9 
and averaged over the intervals 0 — 0.2, 0.2 — 0.4, etc. These 
average quantities were then plotted against 9 and averaged over 
the intervals 0—22.5°, 22.5° — 67-5°, and 67-5° ~ 90° for each 
control ring giving Vn/Vo for each control area. The values of 
v n /Vo for the control areas are given in table III. 

The next step was to determine the slopes of the rectangular 
wedges necessary to cancel the ratio of average finite normal veloc- 
ity to free— stream velocity Vn/v 0 for each control area. The slopes 
(dz/dx)A— 1 , (dz/dx)B— I, and ( ciz/doc) 0 — j for ring I were calculated from 
equations (17), (18) , and (19) . The values of v n /v 0 for the control 
areas of rings II to XVIII were adjusted using equation (20) and the 
factors in table II. Next, the values of (dz/dx)A_jj, (dz/dx)-g_jj, 
and (dz/dx)Q_jj for ring II were determined, and the process was con- 
tinued to ring XVIII, the last ring affecting the wing pressures. 

The slopes of the rectangular wedges are given in table IV. 
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The chord-wise pressure distributions due to the approximate can- 
cellation of the finite velocities normal to the body were finally 
determined using equation (23) . Some consideration was given to 
choosing the chordwise positions at which the pressure coefficients 
of the wing— body juncture were to be calculated. Consider sides A 
and E of the octagonal prism in figure 8 . The rectangular wedges on 
these sides will have leading-edge shock waves which intersect the 
wing surface. For a short distance behind the shock waves a constant 
pressure exists. Calculation shows that the constant pressure region 
covers 6 percent of the chord at the wing-body juncture but only 0.5 
percent at the q/s=1.0 station because of the effect of the tip Mach 
cones of the wedge. At the end of the constant pressure region the 
chordwise pressure distribution has a vertical tangent, and the 
pressure coefficient falls rapidly to a much lover value a short dis- 
tance back. The points at which the pressure coefficients were cal- 
culated were chosen so that they did not coincide with the constant 
pressure peaks at the wing— body juncture. The distance between the • 
points was taken equal to the width of the rings ( 0 . 1 c or 0 . 2 d for 
the present example) to simplify the calculations. 

The coefficients Ci, C z , and C 3 used in equation (23) in cal- 
culating the contributions to the chordwise pressure distributions of 
the first set of rectangular wedges are presented in table V. The 
contributions of the second set of rectangular wedges to the wing 
pressure coefficients were also calculated using equation ( 23 ) and 
table V, taking account of the fact that the second set of wedges 
start a distance 0.2d downstream of the first set. The factors for 
the first set of wedges are Ci=0.4462, C 2 O. 8613 , and Cg=0.4l68 at 
the 0.582 chord position of the 0.05 T}/s station; whereas the 
factors for the second set of rectangular wedges are C 1 = 0 . 6858 , 
02=1.2024, and 03 = 0.5396 at this point. Thus, the contribution of 
all rectangular wedges to the pressure coefficient at any point was 
readily summed. 

The pressure distributions due to canceling approximately the 
finite velocities normal to the body due to the leading— edge sources 
are shown in figure 12 for four of the seven spanwise stations. 

Since the outer stations are farthest from the wedges, the scatter 
of the points decreases and the accuracy improves as the distance 
from the body increases. Curves have been faired through the points 
of the 18— ring solution. For purposes of comparison, the points of 
the chordwise pressure distribution using a 9-ring solution with the 
width of the rings equal to 0.4d have also been plotted. These 
points scatter about the curves for the 18-ring solution, but the 
scatter is not much greater than that of the more accurate method. 
Apparently, the degree of approximation does not change rapidly with 
the number of rings chosen. 
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Over-All Interference Pressure Distributions 
Due to Leading-Edge Sources 

The over-all interference pressure distributions due to the 
leading-edge sources are represented by the sums of the pressure dis- 
tributions of figure 11 plus those of figure 12. The over-all 
pressure distributions for four spanwise stations are shown in figure 
13. This figure brings out a number of significant facts. Eor points 
very near the leading edge of the wing-body juncture, the body is 
effectively an infinite perpendicular reflection plane. Thus, the 
interference pressure coefficient here is 0.52(t/c), just as it is at 
the root chord of the wing alone for either of the leading-edge 
sources. If the body were an infinite perpendicular reflection plane, 
the interference pressure coefficient for all points of the wing-body 
juncture would be 0.52(t/c). However, because the body is curved and 
is not a perfect reflection plane, the pressure coefficient decreases 
rapidly to a much lower value a short distance behind the leading edge 
of the wing-body juncture. 

The over-all interference pressure distributions due to the 
leading— edge sources show increasing positive pressures over the rear 
part of the wing on the inboard sections and over most of the wing on 
the outboard sections because of the influence of the opposite half- 
wing of the combination. The region of influence of one half— wing of 
the wing-body configuration on the body and the other half— wing can be 
determined from the following considerations. Consider the leading 
edge of the wing— body juncture of the left half— wing as a source of 
pressure disturbances which radiate in all directions to the left 
within the Mach cone from the source. A pressure disturbance will 
start to travel in a vertical plane along the intersection of the Mach 
cone with the body at k^° to the wing-body juncture. After this dis- 
turbance has traveled a short distance up the side of the body to a 
new position P, it can be thought of as the source of a number of 
secondary disturbances which radiate in all directions within the 
Mach cone from point P to the left of the plane tangent to the body 
at point P. The" secondary disturbance traveling along the inter- 
section of the Mach cone from P with the body will travel in a 
direction which is k^>° to the element of the cylindrical body through 
P. The forward boundary of the region of influence of the left half- 
wing on the body thus intersects the elements of the cylinder at 45°, 
and would be a straight k$° line if the cylindrical surface of the 
body were unrolled. 

Prom the knowledge of the position of the forward boundary of 
the region of influence of the left half— wing on the body, the 
forward boundary of the region of influence of the left half— wing 
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on the right half— wing may he determined. From the foregoing para- 
graph, it is apparent that the effect of the disturbances starting 
at the leading edge of the left wing— body juncture must travel a 
distance it(d/2) downstream before it can reach the right wing-body 
juncture. For the present example the effect reaches the right 
wing— body juncture at about the 78 - percent chord, and this position 
corresponds to the starting point of the forward boundary of the 
region of influence of the left half— wing on the right half— wing. It 
might be surmised that this boundary is the Mach line on the right 
half— wing which starts at the 78 -percent— chord position of the right 
wing— body juncture. However, this is not the case. Consider secondary 
disturbances from points on the right-hand side of the body. It is 
apparent that some of these secondary disturbances can reach the right 
half— wing along straight paths from their respective sources. When 
account is taken of this effect, it can be shown that the forward 
boundary of the region of influence of the left half— wing on the right 
half— wing actually lies in front of the Mach line from the 7&-percent— 
chord position of the right wing— body juncture. The coordinates of 
the forward boundary are given by the equation 



it + 2 




COS— 1 




2t]+d J J 


( 24 ) 


The | distance is measured positive downstream from the leading 
edge of the wing— body juncture and t] is the lateral distance from 
the wing— body juncture. It is apparent that the forward boundary 
asymptotically approaches a Mach line. 

The forward boundary of the region of influence of the left 
half-wing on the right half— wing, shown in figure 13 , is in close 
agreement with the start of a region of rising positive pressure 
•except near the wing— body juncture. The region of rising positive 
pressure is a consequence of the positive pressure field of the left 
leading— edge source passing around the body. If the present analysis 
had been rigorous to the first order, the region of positive rising 
pressure would probably have started as a cusp at the forward boundary. 
The slight discrepancy noted at the wing-^body juncture is due to the 
approximations of the present method. 

It is apparent from the foregoing discussion of the oven-all 
interference pressure distributions due to the leading-edge sources 
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that no cusps should arise in these pressure distributions in front 
of the forward boundary of the region of influence at the opposite 
half— wing. The cusps in front of this boundary,, shown in figure 13, 
would have been identically canceled if a precise first-order method 
rather than an approximate first— order method had been used to obtain 
the pressure distribution due to cancellation of the finite velocities 
normal to the body. Fob these reasons the cusps in question have 
been -faired as shown by the dashed lines -in figure 13. The effect of 
fairing the first cusps near the Mach line of the leading edge is 
small. However, the fairing of the second and third cusps has an 
appreciable influence on the pressure distributions and was performed 
as follows: The fairing started at the second cusp as a continuation 

of the pressure distribution in front of this cusp and joined the 
pressure distribution on the tangent direction_behind the third cusp. 
The effect of this fairing on the interference pressure drag is con- 
sidered in the next section. 


Interference Pressure Drag of Half— Wings 

As previously pointed out, the interference pressure drag is 
equal to the drag of the wing in combination minus that of the wing 
alone for the configurations considered in this report. The wing 
alone is taken to be the two half— wings joined together. The pres- 
sure drag of the wing alone is considered before determining the 
interference pressure drag. The spanwise distribution of pressure 
drag for the wing alone is shown in figure 14. , This figure is based 
on results of reference 9 and was obtained using the line pressure- 
source theory of R. T. Jones (reference 5). The pressure drag on the 
given half— wing due to its own pressure field accounts for most of 
the wing-alone drag. In fact, for a 10-percent— thick wing the 
pressure drag coefficient is 0.0086, of which 0.0064 is due to the 
pressure fields of the half— wings acting on themselves, while 0.0022 
is due to the pressure fields of the half— wings acting on the opposite 
half— wings. The effect of the wing tips on the spanwise distribution 
of the pressure drag has not been included in the figure, since, as 
previously mentioned, the presence of the tip changes only the distri— - 
bution of the pressure drag within the Mach cones from the leading 
edge of the tips and does not change the over-all wing drag (refer- 
ences 8 and 9) • 

The spanwise variation of the pressure drag of the half— wing 
associated with the over-all interference pressure distribution can 
easily be determined from the over-all interference pressure distri- 
bution due to the leading— edge sources (fig. 13). Denoting the 
average values of P/(t/c) for any spanwise station over the 
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ch.ordv±se interval from —1.0c to — 0.5c as from —0.5c to 0 as 

n 2 , from 0 to 0.5c as hq, and from 0.5c to 1.0c as n 4 , the section • 
drag coefficient due to the over-all interference pressure distribu- 
tion for the leading— edge sources is (n 3 — 114) (t/c) 2 . The 1 section 
drag coefficient for the midchord— line sinks is — 2 (n 2 -H 3 ) (t/c) 2 and 
that for the trailing— edge sources is (ni-n 2 ) (t/c) 2 . Thus the 
resultant section drag coefficient is 

ca = (^ 1 - 3112 + 3113 - 114 ) ( 25 ) 

Equation (25) is valid only for the present wing, "which has a symmet- 
rical double-wedge section with the maximum thickness at the midchord. 

With the aid of equation (25), the spanvise distribution of the 
pressure drag has been determined from the over-all interfer enc e 
pressure distributions of figure 13 , both for the faired and the 
unfaired curves. (See fig. 15 .) Smooth curves have been drawn 
through the calculated points shown in figure 15 as squares for the 
unfaired condition and as circles for the faired condition, 'The 
spanvise distribution of the pressure drag for the two cases is con- 
siderably different. The Importance of the difference can perhaps 
best be assessed by a comparison of the drag coefficients for the 
wing in combination for the two conditions. The drag coefficient of 
either half-wing, neglecting interference effects, is represented by 
the upper curve of figure 14, and, for a 10-percent-thick wing, this 
drag coefficient is 0.0064. The drag coefficient of the wing alone 
is obtained by adding the effect of the opposite half-wing and is 
0.0086 for the same wing thickness. When the drag coefficients 
represented by the two curves of figure 15 are added to 0.0064, the 
resulting drag coefficient for the wing in combination, including 
interference effects, is 0.0075 for the unfaired case and 0.0060 for 
the faired case. Thus, the drag of the wing in combination is 87 
percent of the drag of the wing alone for the over-all interference 
pressure distributions with the cusps, and JO percent without the 
cusps. In either case, the effect of interference is favorable. It 
Is probablp that the 30— percent reduction in drag coefficient repre- 
sents a closer approximation to the first-order answer for the same 
reasons that it is believed the faired over-all Interfer enc e pressure 
distributions represent a closer approximation to the first-order 
pressure distributions than do the unfaired ones. It should be 
noted that the foregoing results apply only to the wing-body con- 
figuration of the example for Mq= and It is conceivable that for 
other configurations a, rid Mach numbers the effect of interference may 
be considerably different. 
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The results considered so far have been for distances of the 
origins of the image sinks from their respective wing— body • junctures 
of 0.05d. The results obtained by the present method should not vary 
much with small changes in the image-sink distances from this value. 
The results of calculations using a 9-^iug scheme and image-sink 
distances of O.ld substantiate this statement. The spanvi.se distri- 
bution of the pressure drag determined from these calculations, using 
the faired curves, are shown in figure 15 as triangles, and the dis- 
tribution is in good accord with that for the 18-ring scheme with 
image— sink distances of 0.053-. 

When the effect on the wing pressures of the pressure field of 
the portion of the body in front of the region of influence of the 
Wing is neglected, the pressure drag of the wing-body configuration 
is less than the sum of the pressure drags of the body alone and the 
wing alone. For a position of the wing well back on the body, the 
effect of the nose should be small; but for a position of the wing 
near the nose, the effect of the nose on the wing may be greater than 
the other interference effects. This appears to be a worthwhile point 
to investigate, especially since it is amenable to simple analysis. 
Another effect which may modify the pressure drag of the wing-body 
combination is the effect that the interference pressure field may 
have on the boundary layer of the combination. 

The analysis has been for Mo= a/ 2. However, this does not 
limit its applicability since any wing-body configuration at a Mach 
number other than */2 has an equivalent wing-body configuration at 
Mq= a/2. To determine the interference pressures on a wing-body 
combination at a Mach number other than a/2, multiply the lo ngitud inal 
dimensions of the wing-body configuration by the factor l/ a/Mq 2 — 1 to 
give an equivalent wing— body combination at Mq= a/ 5. Solve for the 
interference pressures of the equivalent wing— body combination by the 
method of this report. Then the pressure coefficient at a point on 
the actual wing— body configuration at a Mach number other than a/ 5 
is equal to the pressure coefficient at the corresp onding point of 
the equivalent wing-body configuration divided by a/Mq 2 — 1. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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TABLE I 

FACTOR k FOE DETERMINING THE AVERAGE NORMAL 
VELOCITY INDUCED THROUGH THE CONTROL AREAS BY 
THE RECTANGULAR WEDGE COINCIDING WITH SIDE A 
OF THE OCTAGONAL IRISM AND OF WHICH THE REGION 
OF INFLUENCE ON THE BODY BEGINS AT | * 0 


Strip 

Ring 

E 

D 

C 

B 

A 

0 

0 

0 

O.O69 

0.813 * 

I 

0 

0 

0 

.189 

.887 

n 

0 

0 

.062 

.273 

.861 

hi 

0 

.030 

.224 

• 255 

.851 

IT 

.240 

.295 

.218 

.230 

.845 

T 

.242 

.210 

.165 

.223 

.840 

VI 

.180 

.165 

.153 

.218 

.835 

VII 

.160 

.152 

.146 

.214 

.833 

VIII 

.145 

.144 

.142 

.211 

.832 

IX 

.139 

.139 

.139 

.209 

.832 

X 

.136 

.135 

.137 

.208 

.832 

XI 

.134 

.134 

.135 

.207 

.832 

XII 

.133 

.132 

.134 

.206 

.834 

XIII 

.131 

.131 

• 133 

.206 

.836 

XIV 

.130 

.130 

.132 

.206 

.838 

XV 

.130 

.130 

.132 

.206 

.840 

XVI 

.130 

.130 

.132 

.206 

.841 

XVII 

.130 

.130 

• .132 

.206 

.842 

XVIII 
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TABLE H 

FACTORS FOR ADJUSTING THE AVERAGE FINITE NORMAL VELOCITY AT VARIOUS 
CONTROL AREAS DUE TO OCTAGONAL PRISM OF BECTANGULAE WEDGES 
WITH EEGION OF INFLUENCE ON BODY BEGINNING AT | = 0 


-330 

.866 

1.082 


.306 
. 766 
1.015 


Strip 


O.069 ki= 
.813 k z - 
.069 k 3 = 


.285 

1.299 

.285 


• 525 
1.521 
.525 


• 433 
1.412 
.433 


.383 

1.321 

.383 


.386 

1.285 

.366 


.355 

1.261 

.355 


Ring 


=0.813 
= .138 I 
=0 


0.887 
. 378 II 

0 


.851 

.570 IV 
.448 


1.085 

1.050 

.436 


1.082 

.866 VI 
.330 


1.015 
. 7 66 VII 

.306 


• 977 

,710 IS 
,284 


Strip 

Ring 

C 

B 

A 

0.278 
.696 
• 971 

0.348 

1.249 

.348 

0.971 

.696 

.278 

X 

.274 

.686 

.968 

.343 

1.242 

.343 

.968 

.686 

.274 

• 

XI 

.270 
.682 
. 966 

.341 

1.236 

.341 

.966 

.682 

.270 

XII 

.268 

.676 

.967 

.338 

1.235 

.338 

.967 

.676 

.268 

XIII 

.266 

.674 

.967 

-337 

1.233 

.337 

.967 

.674 

.266 

XIV 

.264 

.672 

.968 

.336 

1.232 

.336 

.968 

.672 

.264 

XV 

.264 
.672 
• 970 

.336 

1.234 

.336 

.970 

.672 

.264 

XVI 

.264 

.672 

.971 

.336 

1.235 

.336 

.971 

.672 

.264 

XVII 

.264 

.672 

• 972 

.336 

1.236 

.336 

.972 

.672 

.264 

xvm 
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TABLE III 

RATIOS OF AVERAGE FINITE NORMAL VELOCITY TO FEEE-STREAM VELOCITY 
, v n /V Q THROUGH CONTROL AREAS AFTER CANCELING INFINITE SIDE- 
WASH OF LEADING-EDGE SOURCES FOR 10— HERCENT— THICK WING 


Strip 

Ring 

C 

B 

A 

0 

0 . 

0.001 

I 

0 

.0015 

-.0100 

II 

.0010 

.0065 

-. 0200 

III 

.0080 

.oo4o 

-.0300 

IT 

.0190 

.0070 

-.0350 

V 

.0210 

.0120 

-.0300 

VI 

.0240 

.0110 

-.0350 

VII 

.0275 

.0130 

-.0390 

VIII 

.0300 

.0130 

-.0390 

IX 

.0315 

.0140 

-.o4io 

X 

.0325 

.0140 

-. 0410 

XI 

.0340 

.0140 

-.0420 

XII 

.0340 

.0150 

-.0440 

XIII 

.0350 

.0150 

-.0450 

XIV 

.0360 

.0150 

-.o46o 

XV 

.0370 

.0150 

-.0460 

XVI 

.0380 

.0150 

-.0460 

XVII 

.0390 

.0150 

-. 0470 

XVIII 
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TABLE IV 

SLOPE OF RECTANGULAR WEDGES CANCELING 
AVERAGE FINITE NORMAL VELOCITIES 


(S)„ 

(sx 

(SX 

Ring 

-0.00002 

0.00011 

-0.00125 

I 

.00046 

-.00290 

.01411 

II 

.00131 

-.00917 

.01336 

III 

-.00705 

-.00009 

.01479 

IV 

-.01803 

-.00312 

.00783 

V 

-. 004-99 

-.00453 

-.00601 

VI 

-. 003 U 

.00295 

.01513 

VII 

-.00392 

-.00297 

.00666 

VIII 

•00539 

.00442 

.00063 

IX 

-.00619 

-.ooo 4 o 

.oo 64 i 

X 

-.00707 

-.00585 

-.00721 

XI 

.00292 

.00090 

.00127 

XII 

-.00484 

-.00212 

.00418 

XIII 

.00475 

.00269 

.00344 

XIV 

.00293 

.00445 

.00484 

XV ' 

-.00832 

-.00457 

-.00279 

XVI 

-.00025 

-.00018 

-.00243 

XVII 

-.00468 

-.00176 

-.01358 

XVIII 
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table V 

FACTORS FOR IMHKMIimTG THE WING PRESSURE COEFFICIENTS DUE 
TO THE FIRST SET OF RECTANGULAR WEDGES 


n/s 

= 0 

n/s = 

0.05 

n/s = 

0.125 

Position 

Factor 

Position 

Factor 

Position 

Factor 


Ci = 0 


c x = 0.7788 


Ci = 0.6026 

O.Ol8c 

C 2 * 0 

0.082c 

c 2 = 1.0536 

0.026c 

C 2 - 1.0204 


C3 = 0 


C S = 0 


C3 = 0 


1.0522 


0.4476 


0.3730 

.11 8c 

1.0000 

I .182c 

1.4300 

.126c 

1.3408 

* 

0 


0 


0 


.5324 


.3194 


.2764 

,2l8c 

1.5028 

.282c 

.7780 

.226c 

.6612 


0 


1.0480 


.9824 


.3628 


.2496 


.2212 

.318c 

.9024 

.382c 

1.7356 

.326c 

I.5696 


1.0916 


.8508 


.7372 


.2760 


.68 58 


.6342 

.4l8c ' 

I.8796 

.482c 

1.2024 

.426c 

I.O996 


.9660 


.5396 


.4896 


.7274 


.4462 


.4152 

.518c 

1.2912 

.5820 

,8613 

.526c 

.7980 


.5832 


.4168 


.3840 


.4714 


.3549 


.3322 

.6l8c 

.9144 

.682c 

.6990 

.626c 

.6528 


.4436 


.3443 


.3204 


.3730 


.3000 


.2824 

.718c 

.7368 

.782c 

.5952 

.726c 

.5592 


.3636 


.2553 


.2768 


.3138 


.2616 


.2474 

.8l8c 

.6240 

.882c 

.5212 

.826c 

.4920 


.3100 


.2594 


.2444 


.2730 


.2331 


.2210 

.918c 

.5440 

.982c 

.4648 

.926c 

.4404 


.2712 


.2316 


.2196 


.2422 




.2006 

1.0l8c 

.4832 


_ _ _ 

1.026c 

.4008 


.2412 




.1996 












































































TABLE V.- COBTEJUED 


tj/b - 

O.25 

n/a = 

0.50 

Tj/S - 

0.75 

tj/s = 1.00 

Position 

Factor 

Position 

Factor 

Position 

Factor 

Position 

Factor 


C a = 0.4670* 


0 X - 0.3482 


Cl a 0.2900 


Ci - 0.2536 

-0.065c 

C 2 = .9416 

-0.248c 

C 2 - .8168 

-0.431c 

c 2 - .7280 

-0.6l4o 

C 2 * . 6648 


C 3 = 0 


G3 = 0 


Cg * 0 


& 

11 

0 


0.3036 


0.2346 


0.1981 


0.1746 

.035c 

.9128 

-.l48c 

.6432 

-.331o 

.5284 

-.514c 

.4584 


.3431 


.4076 


.3902 


.3663 


.2322 


.1840 


.1571 


.1394 

.135o 

.5468 

— . 048c 

.4280 

-.231c 

.3639 

— . 4l4c 

.3215 


.8916 


.7648 


.6796 


.6180 


.1898 


.1536 


.1324 


.1182 

•235c 

1.3772 

.052o 

1.1456 

-.131c 

.9678 

-.314c 

.9027 


.6204 


.4928 


.4216 


.37^2 


.5701 


.4852 


.42 99 


.3911 

.335c 

.9750 

• 152c 

.8144 

-.031c 

.7327 

—.21 4c 

.6435 


.4300 


.3564 


.3108 


.2792 


.3759 


.3222 


.2870 


.2611 

.435c 

.7174 

.252c 

.6100 

.0690 

.5455 

~.ll4c 

.4893 


.3434 


.2900 


.2558 


.2312 


.3026 


.2616 

■ 

.2340 


.2136 

.535c 

-5923 

.352c 

.5096 

.169c 

.4566 

— .0i4c 

.4133 


.2898 


.2480 


.2204 


.2004 


.2588 


.2256 


.2026 


.1854 

.625c 

.5114 

.452c 

.4436 

.269c 

.3997 

.086c 

.3634 


.2525 


.2188 


.1958 


.1782 
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TABLE V.- CONCLUDED 


ti/b = 0.25 

*l/s »= 0.50 j 

t]/b = O.75 

r\j s = 1.00 

Position 

Factor 

Position. 

Factor 

Position 

Factor 

Position 

Factor 

0.735c 

0.2280 

.4528 

.2247 

0.552c 

0.2002 

-3964 

.I960 

0.369c 

0.l8o4 

.3579 

.1762 

O.I860 

0,1656 

.3268 

.1614 

.835c 

.2049 

.1*077 

.2027 

.652c 

.1810 

.3592 

.1784 

.469c 

.1639 

.3258 

.1612 

.286c 

.1509 

.2988 

.1481 

• 935c 

.1864 

-3715 

.1852 

.752c 

.1654 

.3296 

.l6i*o 

.5690 

.1505 

.2997 

.1487 

.3860 

.1389 
• 2759 
.1371 

1.035c 

•w<‘ESS» 

mmSEM 

,852c 

.1528 

.3052 

.1520 

.669c 

.1395 

.2783 

.1384 

.486c 

.1291 

.2568 

.1278 

— 

— 

.9520 

.1422 

.2836 

.1416 

.769c 

.1304 

.2602 

.1294 

.5860 

.1210 

.2409 

.1198 

— 

— 

1.052c 

.1332 

.2660 

.1324 

.8690 

.1225 

.2444 

.1217 

.6860 

,1138 

.2268 

.1132 

— 

— 

— 

— 

.969c 

• 1155 
.2308 
.1150 

.7860 

.1076 

.2146 

.1071 

— 

• 

— 

— 

1.069c 

.1094 

.2185 

.1090 

,886c 

.1021 

.2038 

.1016 

— 

— — 

— 

— 

— 


.9860 

.0973 

.1940 

.0968 
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Figure 2. - Sidewash contours for tine pressure source swept 
back 60° Mo = JF. 
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Figure 3. - Source - sink system forming an untapered 
half - wing swept back at 60° with a symmetrical 
double - wedge section. M 0 - . 



38 


MCA EM A9E19 





Figure 4. - Interaction between pressure fields of two half - wings 
separated by half a chord length. M a - 
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Source 


Sink 



I and EE - Leading - edge sources . 
H and JE - Image sources 
M and JZT - Image sinks 


Figure 5.- Source - sink system for reduction of 
infinite sidewash velocities normal to body due 
to leading - edge sources. 
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Figure 6 - Control areas used in approximately canceling finite 
velocities norma! to body. 
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Figure 7, - Source - sink system forming rectangular wedge 
of infinite chord. 



L Contour of octagonal prism 
Body contour 
Plane of the wing 


Figure 8. - Cross - sectional view of regular octagonal prism 
inside of circular body. 
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Figure 9 . - Dimensions of wing - body configuration used 
in illustrative example. 
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Figure 10. - Pressure distribution at right wing- body juncture due to 
leading - edge source of left half - wing and due to reduction of 
infinite side wash vebcities of leading - edge sources. Mo=JF. 







Figure //. - Spanwise variation of pressure distribution due to leading - edge source of left half - wing 
and due to reduction of infinite s/dewash velocities of leading -edge sources. 




(a) y/S - 0 





(b) v /s -- .125 


Figure 12 interference pressure distributions due to approximate cancellation 
of finite velocities normal to body due to leading - edge sources. 


S 


6TSC6V MS TO VCt 




Percent chord 
(c) v /$ - .50 % 



(d) v /s = LOO 


Figure 12. - Concluded. 
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6ia6V HH VOW 





£ 


, 1?/s ‘ LOO 



Figure /j,- Over -all interference pressure distributions due to leading -edge sources for wing in combination 
with body. 
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Lateral distance , fraction of semi span , y/s 


Figure 14 Spanwise distribution of the pressure drag for either half 
wings joined at root. M n =[2^. 
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Figure 15.- Spanwise distribution of the pressure drag associated with the 
over- all interference pressure distribution. 
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